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(— I , Abstract 

I ' The anomalous effective action describing the coupUng of gravity to a non-abehan gauge theory 

flj ' can be determined by a variational solution of the anomaly equation, as shown by Riegert long 

1-^ . ago. It is given by a nonlocal expression, with the nonlocal interaction determined by the Green's 

function of a conformally covariant operator of fourth order. In recent works it has been shown 
\^ I that this interaction is mediated by a simple pole in an expansion around a Minkowski background, 

rf-N ' coupled in the infrared in the massless fermion limit. This result relies on the local formulation of 

C^^ ' the original action in terms of two auxiliary fields, one physical scalar and one ghost, which take the 

role of massless composite degrees of freedom. In the gravity case, the two scalars have provided 
jy~, , ground in favour of some recent proposals of an infrared approach to the solution of the dark energy 

<^^ ' problem, entirely based on the behaviour of the vacuum energy at the QCD phase transition. As 

a test of this general result, we perform a complete one- loop computation of the effective action 
describing the coupling of a non-abelian gauge theory to gravity. We confirm the appearance of an 
anomaly pole which contributes to the trace part of the TJJ correlator and of extra poles in its 
trace-free part, in the quark and gluon sectors, describing the coupling of the energy momentum 
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C^ ' tensor (T) to two non abelian gauge currents (J). 
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1 Introduction 

The study of the effective action describing the coupUng of a gauge theory to gravity via the trace 
anomaly [T] is an important aspect of quantum field theory, which is not deprived also of direct 
phenomenological implications. This coupling is mediated by the correlator involving the energy 
momentum tensor together with two vector currents (or TJJ vertex), which describes the interaction 
of a graviton with two photons or two gluons in QED and QCD, respectively. At the same time, 
the vertex has been at the center of an interesting case study of the renormalization properties of 
composite operators in Yang Mills theories p], in the context of an explicit check of the violation of 
the Joglekar-Lee theorem [3] on the vanishing of S-matrix elements of BRST exact operators. In this 
second case it was computed on-shell, but only at zero momentum transfer. In this work we are going 
to extend this computation and investigate the presence of massless singularities in its expression. 
These contribute to the trace anomaly and play a leading role in fixing the structure of the effective 
action that couples QCD to gravity. The analysis of [2], which predates our study, unfortunately does 
not resolve the issue about the presence or the absence of the anomaly pole in the anomalous effective 
action of QCD because of the restricted kinematics involved in that analysis of the TJJ vertex, and 
for this reason we have to proceed with a complete re-computation. 

Anomaly poles characterize quite universally (gravitational and chiral) anomalous effective actions, 
in the sense that account for their anomalies. They have been identified and discussed in the abelian 
case both by a dispersive analysis [4] and by a direct explicit computation of the related anomalous 
Feynman amplitudes quite recently [5l [6] . Extensive analysis in the case of chiral gauge theory for 
anomalous f/(l) models have shown the close parallel between solutions of the Ward identities, the 
coupling of the poles in the ultraviolet and in the infrared region and the gravity case [T, "5]. 

It is therefore important to check whether similar contributions appear also in non-abelian gauge 
theories coupled to gravity. We recall that the same pole structure is found in the variational solution 
of the expression of the trace anomaly, where one tries to identify an action whose energy momentum 
tensor reproduces the trace anomaly. This action, found by Riegert long ago [S], is nonlocal and 
involves the Green function of a quartic (conformally covariant) operator. The action describes the 
structure of the singularities of anomalous correlators with any number of insertions of the energy 
momentum tensor and two photons (T^JJ), which are expected to correspond both to single and to 
higher order poles, for a sufficiently high n. For obvious reasons, explicit checks of this effective action 
using perturbation theory - as the number of external graviton lines grows - becomes increasingly 
difficult to handle. The TJJ correlator is the first (leading) contribution to this infinite sum of 
correlators in which the anomalous gravitational effective action is expanded. 

Given the presence of a quartic operator in Riegert 's nonlocal action, the proof that this action 
contains a single pole to lowest order (in the TJJ vertex), once expanded around fiat space, has been 
given in [^ by Giannotti and Mottola, and provides the basis for the discussion of the anomalous 
effective action in terms of massless auxiliary fields contained in their work. The auxiliary fields 



are introduced in order to rewrite the action in a local form. We show by an explicit computation 
of the lowest order vertex that Riegert's action is indeed consistent in the non-abelian case as well, 
since its pole structure is recovered in perturbation theory, similarly to the abelian case. Therefore, 
one can reasonably conjecture the presence of anomaly poles in each gauge invariant subsets of the 
diagrammatic expansion, as the computation for the non-abelian TJJ shows (here for the case of 
the single pole). In particular, this is in agreement with the result of \^, where it is shown that, 
after expanding around flat spacetime, the quartic operator in Riegert's action becomes a simple l/D 
nonlocal interaction (for the TJJ contribution), i.e. a pole term. We remark that the identification 
of a pole term in this and in others similar correlators, as we are going to emphasize in the following 
sections (at least in the case of QED and for the sector of QCD mediated by quark loops), requires 
an extrapolation to the massless fermion limit, and for this reason its interpretation as a long-range 
dynamical effect in the gravitational effective action requires some caution. In QCD, however, there 
is an extra sector that contributes to the same correlator, entirely due to virtual loops of gluons in 
the anomaly graphs, which remains unaffected by the massless fermion limit. The appearance of such 
a singularity in the effective action, however, does not necessarily imply that its contribution survives 
in the physical S-matrix. We will also establish the appearance of other singularities in the trace-free 
form factors which, obviously, are not part of Riegert's action. 

We will comment in our conclusions on the possible implications of these results and on some 
recent proposals to link this type of behaviour [lOl [11] to cosmology and to the dark energy problem. 
We also remark that, in general, the coefficient in front of the trace anomaly, for a given theory, can be 
computed in terms of its massless fields content, and as such it is well known. However, the structure 
of the effective action and the characterization of its fundamental form factors at nonzero momentum 
transfer and its complete analytical structure is a novel result. In this respect, the classification of all 
the relevant tensor structures which appear in the computation of this correlator is rather involved 
and has been performed in the completely off-shell case. We remark that the complexity of the final 
expression, in the off-shell case, prevents us from presenting its form. For this reason we will give only 
the on-shell version of the complete vertex, which is expressed, as we have mentioned, only in terms 
of three fundamental form factors. 

Concerning the phenomenological relevance of this vertex, we just mention that it plays an essential 
role in the study of NLO corrections to processes involving a graviton exchange. In fact, in theories with 
extra dimension, where a low-gravity scale and the presence of Kaluza-Klein excitations may enhance 
the rates for processes mediated by gluons and gravitons, the vertex appears in the hard scattering 
of the corresponding factorization formula [12j and has been computed in dimensional regularization. 
However, to our knowledge, in all cases, there has been no separate discussion of the general structure 
of the vertex (i.e. as an amplitude) nor of its Ward identities, which, in principle, would require a more 
careful investigation because of the trace anomaly. Anomalous amplitudes, in fact, are defined by the 
fundamental Ward identities imposed on them, that we are going to derive from general principles. We 
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Figure 1: The diagrams describing the anomaly pole in the dispersive approach. Fig. (a) depicts the singularity 
of the spectral density p{s) as a spacetime process. Fig. (b) describes the anomalous pole part of the interaction 
via the exchange of a pole. 



cover this gap and show, that both dimensional regularization and dimensional reduction reproduce 
the correct Ward identity satisfied by this vertex, showing at the same time that the use of these 
regularizations is indeed appropriate. Results for this vertex will be given only in the on-shell case, 
since in this case the result can be expressed in terms of just three form factors. We have computed 
also the off-shell effective action, but its expression is rather lengthy and will not be discussed here, 
since it is gauge dependent and of less significance compared to the on-shell result. Most of our work 
is concerned with a technical derivation of the leading contribution to the anomalous effective action 
of QCD coupled to gravity. We have summarized in our conclusions a brief discussion of the relevance 
of this study in the ongoing attempt to link the trace anomaly and QCD to a possible alternative 
solution of the problem of dark energy, using this effective action as an intermediate step |lUt [TT] . 

2 Anomalous effective actions and their variational solutions 

In this section we briefly review the topic of the variational solutions of anomalous effective actions, 
and on the local formulations of these using auxiliary fields. 

One well known result of quantum gravity is that the effective action of the trace anomaly is given 
by a nonlocal form when expressed in terms of the spacetime metric g^u ■ This was obtained [9] from 
a variational solution of the equation for the trace anomaly [I] 

Tli = bF + b' (E-^nRj +b"UR + cF''^"'F^^, (1) 

(see also |131 ^^ for an analysis of the gravitational sector) which takes in I? = 4 spacetime dimensions 
the form 
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Here, the parameters b and b' are the coefficients of the Weyl tensor squared, 

F = Cx^.pC^'"''' = Rx^upR^^"'" - 2R^,R^^ + ^ (3) 

and the Euler density 

E =*Rxpup "R^^"' = Rxp.pR''^'"' - ^Rp.R"" + i?" (4) 

respectively of the trace anomaly in a general background curved spacetime. Notice that the last term 
in d^l) is the contribution generated in the presence of a background gauge field, with coefficient c. 
For a Dirac fermion in a classical gravitational {gfj.u) and abelian (Aa) background, the values of the 
coefficients are b = 1/(320 vr^), and b' = — ll/(57607r^), and c = — e^/(247r^), with e being the electric 
charge of the fermion. One crucial feature of this solution is its origin, which is purely variational. 
Obtained by Riegert long ago, the action was derived by solving the variational equation satisfied 
by the trace of the energy momentum tensor. Aj (x, x') denotes the Green's function inverse of the 
conformally covariant differential operator of fourth order, defined by 

V„, (VV + 2R^'•' - iRgf^A V^ = D^ + 2i2^^V/,V^ + ^(V^'i?)V^ - ^RO . (5) 



Given a solution of a variational equation, it is mandatory to check whether the solution is indeed 
justified by a perturbative computation. One specific feature of these solutions is the presence of 
anomaly poles. In previous works we have elaborated on the significance of these interactions, extracted 
from a direct perturbative computation, by a painstaking analysis of anomaly graphs under general 
kinematical conditions, and not just by a dispersive approach. The dispersive approach allows to 
connect this behaviour of the spectral density to a very specific infrared configuration. 

2.1 The kinematics of an anomaly pole 

In our conventions we will denote with p and q the outgoing momenta of the two photons/gluons 
and with k the incoming momentum of the graviton. s = {p + qf' denotes the invariant mass of the 
external graviton line. A computation of the spectral density p{s) of the TJJ amplitude in QED 
shows that this takes the form p{s) ~ (^(s). The configuration responsible for the appearance of a 
pole is illustrated in Fig. [1] (a). It describes the decay of a graviton line into two on-shell photons. 
The decay is mediated by a collinear and on-shell fermion-antifermion pair and can be interpreted as 
a spacetime process. The corresponding interaction vertex, described as the exchange of a pole, is 
instead shown in Fig. [T] (b). The actual process depicted in Fig. [T] (a) is obtained at diagrammatic 
level by setting on-shell the fermion/antifermion pair attached to the graviton line. This configuration, 
present in the spectral density of the diagram only for on-shell photons, generates a pole contribution 
which can be shown to be coupled in the infrared. This means that if we compute the residue of the 
amplitude for s — )• we find that it is non- vanishing. In the general expression of the vertex, a similar 



configuration is extracted in the high energy limit, not by a dispersive analysis, but by an explicit 
(off-shell) computation of the diagrams. Clearly, the pole, in this second case, has a vanishing residue 
as s — )• 0, but is nevertheless a signature of the anomaly at high energy. Either for virtual or for real 
photons, a direct computation of the vertex allows to extract the pole term, without having to rely 
on a dispersive analysis. This point has been illustrated in our previous computations of the chiral 
anomaly vertex [5] and in the computation of the T J J vertex for QED [S] . The identification of this 
singularity in the case of QCD is in perfect agreement with those previous results. 



2.2 The single pole from A4 

In the case of the gravitational effective action, the appearance of the inverse of A4 operator seems to 
be hard to reconcile with the simpler l/D interaction which is predicted by the perturbative analysis 
of the TJJ correlator, which manifests a single anomaly pole. In |^, Giannotti and Mottola show step 
by step how a single pole emerges from this quartic operator, by using the auxiliary field formulation 
of the same effective action. Clearly, more computations are needed in order to show that the nonlocal 
effective action consistently does justice of all the poles (of second order and higher) which should 
be present in the perturbative expansion. Obviously, the perturbative computations - being either 
based on dispersion theory or on complete evaluations of the vertices, as in our case - become rather 
hard as we increase the number of external lines of the corresponding perturbative correlator. For 
instance, this check becomes almost impossible for correlators of the form TTTT or higher, due to the 
appearance of a very large number of tensor structure in the reduction to scalar form of the tensor 
Feynman integrals. In the case of TJJ the computation is still manageable, since it does not require 
Feynman integrals beyond rank-4. 

Expanding around flat space, the local formulation of Riegert's action, as shown in [UlIS], can be 
rewritten in the form 

which is valid to first order in the fluctuation of the metric around a flat background, denoted as h^u 

(7) 



9^J.l' = ViJiu + k/i, 



^lv^ 



VIGttGat 



with Gn being the 4-dimensional Newton's constant. The formulation in terms of auxiliary fields of 
this axion gives 
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.^'U^-^i,' + ^F^pF^P^ 



(8) 



where (p and ifj are the auxihary scalar fields. They satisfy the equations 

V^'^^DV, (9) 

^v^' = |F,^F-^ (10) 

n<^ = -|. (11) 

In order to make contact with the TJJ amplitude, one needs the expression of the energy momen- 
tum extracted from ([8]) to leading order in /i^,y, or, equivalently, from ([6]) that can be shown to take 
the form 

Ta'^omi^) = I ia^-'a - d^dn^ J d^x' D-i, [f^^f-^] ^^ . (12) 

Notice that Tanom is the expression of the energy momentum tensor of the theory in the background 
of the gravitational and gauge fields. We recall, in fact, that in the QED case, for instance, the energy 
momentum tensor of the theory is split into the free fermionic part Tj, the interacting fermion-photon 
part Tfp and the photon contribution Tp/^ which are given by 

T^" = -i^^(^''^''^ + g>"'{i'il^j^dx4^ - WV'), (13) 



and 



where the current is defined as 



T!f^ = - eJ^^'A"^ + eg'^'J^Ax , (14) 

T^h = F^^F'x - \9'"'F^''Fxp, (15) 



JM(2;) = V'(x)7^V(2;) . (16) 



The connected components of TJJ can be obtained directly from the quantum average of Tp, 
defined as the sum of the fermion contribution and its interaction part with the photon field, 

T^^ = Tf + T'}:^. (17) 

In the formalism of the background fields, the TJJ correlator then can be extracted from the 
defining functional integral 

via two functional derivatives respect to the background field A^ and generates the effective action 



^'''''y>-6A^{x)SAp{y) 



^ anom '^ ^ ' \^^J 
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We have separated in (fT9|) the pole contribution Tanom from the rest of the amphtude (F), which 
does not contribute to the trace part. Notice that Tanom, derived from either the classical generating 
functional (J12p given by Riegert's action or from the direct perturbative expansion of (J19p . should 
nevertheless coincide, for the pole term not to be a spurious artifact of the variational solution. In 
particular, a computation performed in QED shows that the pole term extracted from Tanom via 
functional differentiation 

Knom{P,Q) - Jdxjdye SA^{x)A^iy) ' IStt^ k^ ^^ " '^^J^ l^'?) 

(20) 

with 

n°^(p,g)^(p -5)5"^ -<?-/, (21) 

indeed coincides with the result of the perturbative expansion, as defined from the first term on the 
rhs of (|19|) . Thus, the entire contribution to the anomaly is extracted form Tanom as 

As we have already mentioned, the full action ([2]), varied several times with respect to the background 
metric g^u and/or the background gauge fields Aq. gives those parts of the correlators of higher order, 
such as {TTT...JJ) and (TTT...), which contribute to the trace anomaly. In particular, the anomalous 
contributions of the T" J J's vertices are obtained by varying the local action both respect to the metric 
and to the gauge fields. 

3 The energy momentum tensor and the Ward identities 

Moving to the QCD case, we introduce the definition of the QCD energy-momentum tensor, which is 
given by 

T^u = -5^.£QCD-F;pF7-^9^.5'K5''^^) + |K5,x5'^^^) + A^a,(a'^yi^)) 



g^uT^nom = cF^pF-P = -^F^pF-P. (22) 
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V'7^(^. - igT'^ADi; - ^pCb, + igT'^A^j^i; + V7.(^m " ^5^"^")^' 



^{d^ + igT''Al)jM + d^iJ'^id^u;'' - gP^'AlJ) + d.uf'id^u^ - g P'"' Alco'') , (23) 



where F!^^^ is the non-abelian field strength of the gauge field A 

f;, = d,A: - d^Al + gr'^'A'^Al (24) 



and we have denoted with a;° the Faddeev-Popov ghosts and with lo"" the antighosts, while ^ is the 
gauge-fixing parameter. The T"'s are the gauge group generators in the fermion representation and 



f are the antisymmetric structure constants. For later use, it is convenient to isolate the gauge-fixing 
and ghost dependent contributions from the entire tensor 



T^u^- = -^[A'^Md-A-) + A^^d,{d-A-)]--^g,, 



(25) 
r^t = d^oj'^Dfu' + d^oJ'^Dfu'-g^.d^O^Dfu;'. (26) 



--{d-Af + dfiA^^d-A^) 



The coupling of QCD to gravity in the weak gravitational field limit is given by the interaction 
Lagrangian 



Cint = --^h''''T^,. (27) 
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2 

Notice that T^jy as defined in Eq. (I23p is symmetric, while traceless for a massless theory. The 
symmetric expression can be easily found as suggested in [16j, by coupling the theory to gravity and 
then defining it via a functional derivative with respect to the metric, recovering (j23|) in the flat 
spacetime case. 

The conservation equation of the energy-momentum tensor takes the following form off-shell |17tll8j 

SS^ , ^ -5S K„ /55 ^ . SS\ ^ ^^ SS 



n- - T — d^oj - OvOJ - — 



+ dA^iTl^]-—duU^''-duU^''—, (28) 



where a^i, = ^[7^, 7j/]. It is indeed conserved by using the equations of motion of the ghost, antighost 
and fermion/antifermion fields. The off-shell relation is particularly useful, since it can be inserted 
into the functional integral in order to derive some of the Ward identities satisfied by the correlator. 
In fact, the implications of the conservation of the energy-momentum tensor on the Green's functions 
can be exploited through the generating functional, obviously defined as 



Z[J,r],fj,x,X,h] = / VAVil)V'ipVojVu) expli d!^x{C + J^A^" 



.T'^")), 



+^^ + Vry + xw + cjx + V^"") | > (29) 

where C is the standard QCD action and we have added the coupling of the energy-momentum tensor 
of the theory to the background gravitational field h^^, which is the typical expression needed in the 
study of QCD coupled to gravity with a linear deviation from the flat metric. We have denoted with 
J, 77, 77, X, X the sources of the gauge field A (J), the source of the fermion and antifermion fields (ry, rf) 
and of the ghost and antighost fields (x, x) respectively. The generating functional W of the connected 
Green's functions is, as usual, denoted by 

exp ^ W[J, ry, r?, x, X, h] = ^[qooooq] (30) 
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(normalized to the vacuum functional) and the effective action, defined as the generating functional 
r of the 1-particle irreducible and truncated amplitudes. This is obviously obtained from 14/^ by a 
Legendre transformation respect to all the sources, except, in our case, h^^, which is taken as a 
background external field 

r[Ac, V^c, V'c, ^c, Wc, h] = W[J, r], f], X, X,h]- / d'^x {J^A^ + t/V'c + ipcV + X^c + uicx) ■ (31) 

The source fields are eliminated from the right hand side of Eq. (|3ip inverting the relations 

,, 6W , 5W J 6W 6W _ 6W 

oJ/i or] dr] dx ox 

so that the functional derivatives of the effective action T with respect to its independent variables are 

<5r <5r _ dV 6T _ 6T 

--^t^^ TT = -V, T^ = -V, -^ = -X, T^ = -X, (33) 
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Oipc 
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■ -X, 


<5r 

6Qc 


^ -X' 


sr 5w 










5h^y 5h^y 











and for the source /i^jy we have instead 

XT Xf^/ 

(34) 

The conservation of the energy-momentum tensor summarized in Eq. ()28p in terms of classical fields, 
can be re-expressed in a functional form by a differentiation of W with respect to h^i, and the use of 
Eq. (j28p under the functional integral. We obtain 



+ xdu-^ + du^X, (35) 

Sx ^X 

and finally, for the one particle irreducible generating functional, this gives 

- ^d^u^.-d'^ojc^, (36) 

OUJr OUJc 



obtained from Eq. (j35|) with the help of Eqs. (|32l (|33j) . (|M|) . We summarize below the relevant Ward 
identities that can be used in order to fix the expression of the correlator. 

• Single derivative general Ward identity 

The Ward identities describing the conservation of the energy-momentum tensor for the one- 
particle irreducible Green's functions with an insertion of T^i^ can be obtained from the functional 
equation (j36p by taking functional derivatives with respect to the classical fields. For example. 
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the Ward identity for the graviton - gluon gluon vertex is obtained by differentiating Eq. (|36p 
with respect to A^^{xi) and A'^a{x2) and then setting ah the external fields to zero 

+ d^ [gav^'^ixi - x)D~^{x2,x) + gj3v5'^{x2 - x)D~l^{xi,xyj 

(37) 

where D~a{xi,X2) is the inverse gluon propagator defined as 

D-l{xi,X2) = {Aa,{xi)Ap{x2))trunc = ^ ^ , ^^ a (38) 

M°(xi)M?(x2) 

and where we have omitted, for simplicity, both the colour indices and the symbol of the T- 
product. The first Ward identity (j37p becomes 

k^ {Tf,u{k)Ao,{p)Ap{q))trunc = qfiDa^{p)gi5iy + PfJ,D^^{Q)gau - QiyDa^ip) " PuD^lio) ■ (39) 

Trace Ward identity at zero momentum transfer 

It is possible to extract a Ward identity for the trace of the energy-momentum tensor for the 
same correlation function using just Eq. (I39p . In fact, differentiating it with respect to p^ (or 
q^) and then evaluating the resulting expression at zero momentum transfer {p = —q) we obtain 
the Ward identity in d spacetime dimensions 

{Tli{^)A^{p)Ap{-p))trunc={2-d + p-^^D-l{p) (40) 

where the number 2 counts the number of external gluon lines. For d = 4 and using the 
transversality of the one-particle irreducible self-energy, namely 

^a/3 W = (/5a/3 " Paq^mp"), (41) 

the Ward identity in Eq. ()40p simplifies to 



{Tl^{Qt)Ao.{p)Ap{-p))trunc = 2p\p^ga^ - Paqp)-p^{p'). (42) 



dli 

' dp^ 

The trace Ward identity in Eq. (j40p at zero momentum transfer can also be explicitly related to 
the /3 function and the anomalous dimensions of the renormalized theory. These enter through 
the renormalization group equation for the two-point function of the gluon. Defining the beta 
function and the anomalous dimensions as 

OQ O I OTfl 

^(5) = P-Q-^ 7(5) = M^log a/Z^, m-frn{g) = l^-Q- (43) 

11 



and denoting with Z^ the wave function renormahzation constant of the gluon field, with g 
the renormalized couphng, and with m the renormahzed mass, the trace Ward identity can be 
related to these functions by the relation 



(:rm)Ao.(ji)Afi[-v))trunc 



/?(5)^-27(5)+m(7„^(5)-l)^ 



^->\ 



(44) 



• Two-derivatives Ward identity via BRST symmetry 

We can exploit the BRST symmetry of the gauge-fixed lagrangian in order to derive some general- 
ized Ward (Slavnov- Taylor) identities. We start by computing the BRST variation of the energy- 
momentum tensor, which is given by 






-igtpfXuj"-, 



(45) 
(46) 

(47) 

(48) 
(49) 



where A is an infinitesimal Grassmann parameter. 

A careful analysis of the energy-momentum tensor presented in Eq. (I23p shows that the fermionic and 
the gauge part are gauge invariant and therefore invariant also under BRST. Instead the gauge-fixing 
and the ghost contributions must be studied in more detail. Using the transformation equations (j45p 
and dm in ([2S|) one can prove the two identities 



A T^--^- 



-Ald^SCj'^ - A^duSu^ + g, 



]-d ■ A^Soj" + A^pdPSoj" 



(50) 
(51) 



which show that the ghost and the gauge-fixing parts of the energy-momentum tensor (times the 
anticommuting factor A) can be written as an appropriate BRST variation of ghost/antighost and 
gauge contributions. Their sum, instead, can be expressed as the BRST variation of a certain operator 
plus an extra term which vanishes when using the ghost equations of motion 



A (t^V- + 7]?') = ^ -0,0'^ a: - d^uj'^Al + g,, I^A'^pdpO'' + \d ■ A'^a;") 



+ g^,\xiI:''dPDfuj\ (52) 



which shows explicitly the structure of the gauge- variant terms in the energy-momentum tensor. Using 
the nilpotency of the BRST operator (5^ = 0), the BRST variation of T^y is given by 



A 



5T,, = 6{T^J- + T^!:) = - [A';,d,df>Dfco' + Ald^d^D 



ah b 
P ^ 



g^^d'^iAyPDfu') 



(53) 
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where it is straightforward to recognize the equation of motion of the ghost field on its right-hand 
side. Using this last relation, we are able to derive some constraints on the Green's functions involving 
insertions of the energy-momentum tensor. In particular, we are interested in some identities satisfied 
by the correlator (Tf^^A'^A'a) in order to define it unambiguously. For this purpose, it is convenient 
to choose an appropriate Green's function, in our case this is given by {T^yd°^ A'^ufi) , and then exploit 
its BRST invariance to obtain 

5{T^,d''Alu'>) = {ST^.d'^Aloj'') + A(r^,a"Z)>^L^^) - hr^.d^'Ald^ A''^) = 0, (54) 

where the first two correlators, built with operators proportional to the equations of motion, contribute 
only with disconnected amplitudes, that are not part of the one-particle irreducible vertex function. 
From Eq. (|54|) we obtain the identity 

9^lC(^M-(^)^"(^l)4(^2))t™nc = 0, (55) 

which in momentum space becomes 

p"/(r^,(A:)^^(p)^^(g))t™„, = 0. (56) 

A subtlety in these types of derivations concerns the role played by the commutators, which are 
generated because of the T-product and can be ignored only if they vanish. In general, in fact, the 
derivatives are naively taken out of the correlator, in order to arrive at Eq. (j56|) and this can generate 
an error. In this case, due to the presence of an energy momentum tensor, the evaluation of these 
terms is rather involved. For this reason one needs to perform an explicit check of (j56p to ensure the 
consistency of the formal result in a suitable regularization scheme. As we are going to show in the 
next sections, these three Ward identities turn out to be satisfied in dimensional regularization. 

4 The perturbative expansion 

The perturbative expansion is obtained by taking into account all the diagrams depicted in Figs. O 
[3l m where an incoming graviton appears in the initial state and two gluons with momenta p and q 
characterize the final state. The different contributions to the total amplitude are identified by the 
nature of the internal lines and are computed with the aid of the Feynman rules defined in Appendix 
lAl Each amplitude is denoted by F, with a superscript in square brackets indicating the figure of the 
corresponding diagram. 

The contributions with a massive fermion running in the loop are depicted in Fig.[2l for the triangle 
in Fig. [2^ we obtain 

- 'I rS'(P. ,) = -1/ ^iT^rn I i^t. {km - ,. ; + p)^_i_,,^-J-,„^-^} (57) 
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exch. 



(d) 



Figure 2: The fermionic contributions with a graviton /i^^ in the initial state and two gluons A'^,A'p in the 
final state. 



A" h 




A» 



(a) 





Figure 3: The gauge contributions with a graviton h^i, in the initial state and two gluons A'^,A% in the final 



state. 



where the color factor is given by Tr(T^T") = ^5"''; the diagram in Fig. ^ contributes as 



if rS'(p.<rt = -f»'^(^"^')/ (£f''{^^-7-^ 



-7/3- 



ip 1 f 1 (58) 

m ji — m ] 

with the vertices V"^j^(/ — q,l+p) and W'^^^ defined in Appendix|Al Eqs. (jl06p and (|109p respectively. 
The remaining diagrams in Fig. [2] are obtained by exchanging a •(-)• /3 and p -^ q 



^f r|J'j^''(p,g) 



-^f rjfj2^*(p,^) 



-^ f rjfj'i?^^' 9) 



-^f rjf;i^^(p,g) 



a <->• /3 
p^ q 



p^ q 



(59) 
(60) 



Moving to the gauge sector we find the four contributions in Fig. O the first one with a triangular 
topology is given by 

dH 1 



'^Ir%if{p,q) = -l9^r'^^f'^l^ 



Yl^{l+pf{l-qY 



VG99{l-q,-l-p) X 



V^„^{-Ul+p,-p) V^f^{-l + q,l, 



(61) 
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where the color factor is f^'^^f"'^^ = Ca u . Those in Figs. [8)3 and [3}:, containing gluon loops attached 
to the graviton vertex, are called "t-bubbles" and can be obtained one from the other by the exchange 
of a o /? and p ■h^ q. The first "t-bubble" is given by 

_ . ^ ^[3b] a6 . . _ _ 1 /^ 2 .ade Me f dH VJ^jfl^H, I - P, -<l) VplA^, -p, -I + p) 

which is multiplied by an additional symmetry factor 2 • There is another similar contribution obtained 
from the previous one after exchanging a -f-)- /3 and p ■^ q 



^2 V^"/3 (^'«) = -^2 I'^Lp (P^l) 



a -n- /3 
p^ q 



(63) 



The last diagram with gluons running in the loop is the one in Fig. [3li which is given by 



^2 V«/3 (^'9) =22^7 (2^ lHl-p-1? ' ^''^ 



where V^ is the four gluon vertex defined as 



/lupa 



r'^f^^gppg.a - g^agup) + r''f'''{g^ugp. - g^agup) 

+ r'''f''(.gp.gpa-g,pg.a)] (65) 

and therefore 

^'^/y^J = -CAS''%\,^ = -CA6''''{ga.gpp + gapg^a-2g^0.p), (66) 

so that the amplitude in Eq. (j64|) becomes 



4 rLir (p. «) = --s'c.s'- f ^ ''''^"''-'•'-^-f^'^ . (67) 



-2 V-.,5 ^'■•'" - -52" "'-»" J (2^ W^ 

In the expression above we have explicitly isolated the color factor C^5"^ and the symmetry factor ^. 
Finally, the ghost contributions shown in Fig. U] are given by the sum of 



'2 V^«/3 Vi^'^^ 2^ J ^ y (27r)4 p(/+p)2(/_g) 

for the triangle diagram in Fig. [J^ and 

_ ;,'^ T^"^^^"-^ ( rn n\ — '^ n'2 fade fbde f ^ I C pyggr {I - q)i3 

^2^>^-^P^P^'i>- 2^ ^ / y (27r)4 f{l-q)^ ^^^> 

for the "T-bubble" diagram shown in Fig. dJ:. The two exchanged diagrams are obtained from those 
in Eqs. ([68]) and ([69|) with the usual replacement 0/3 and p ^ q. 



' 2 Ml'"/? ^P'^' ~ 2 ^'^"'^ 



* 2 Ml'"/? '^^' ^'' ~ 2 ^'^"'^ 
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Figure 4: The ghost contributions with a graviton /i^t/ in the initial state and two gluons A°,ylo in the final 



state. 



Having identified the different sectors we obtain the total amplitude for quarks, denoted by a "q" 
subscript 

r^!',.„,(p, g) = v^llt^v. 9) + rlS;^p, q) + rif;] ;^p, ,) + r[fj;;^p, ,) (72) 

and the one for gluons and ghosts as 

rt,ua,iP, ^) = E [rJil?(P' 1) + rS'l?(p, 9) + tJI^Cp, g) + rj,l;^p, g)] . (73) 

i=3,4 

5 The on-shell (TAA) correlator, pole terms and form factors 

We proceed with a classification of all the diagrams contributing to the on-shell vertex, starting from 
the gauge invariant subset of diagrams that involve fermion loops and then moving to the second set, 
the one relative to gluons and ghosts. The analysis follows rather closely the method presented in 
the case of QED in previous works [H [5], with a classification of all the relevant tensor structures 
which can be generated using the 43 monomials built out of the 2 of the 3 external momenta of the 
triangle diagram and the metric tensor g^i,. In general, one can proceed with the identification of a 
subset of these tensor structure which allow to formulate the final expression in a manageable form. 
The fermionic triangle diagrams, which define one of the two gauge invariant subsets of the entire 
correlator, can be given in a simplified form also for off mass-shell external momenta, in terms of 
13 form factors as in [4, 5J while the structure of the gluon contributions are more involved. Some 
drastic semplifications take place in the on-shell case, where only 3 form factors - both in the quark 
and fermion sectors - are necessary to describe the final result. 
We write the whole amplitude T'^'^"" (p, q) as 

r'^-/3(p, q) = rp/5(p, q) + rf -/'(p, g), (74) 

referring respectively to the contributions with quarks (Fg) and with gluons/ghosts (F^) in Eqs. ()72p 
and (|73|) . We have omitted the color indices for simplicity. The amplitude F is expressed in terms of 
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3 tensor structures and 3 form factors renormalized in the MS scheme 



(75) 



i=l 



One comment concerning the choice of this basis is in order. The 3 form factors are more easily 
identified in the fermion sector after performing the on-shell hmit of the off-sheh amphtude, where the 
13 form factors introduced in [H |5] for QED simphfy into the 3 tensor structures that will be given 
below. It is then observed that the tensor structure of the gluon sector, originally expressed in terms 
of the 43 monomials of [H [5], can be arranged consistently in terms of these 3 reduced structures. 
The tensor basis on which we expand the on-shell vertex is given by 






-2u''^ip,q)[sg'''' + 2ip''p'' + q''q'')-4{p''q'' + q^'p'')], 
{p^q^+p^q^)9' 



(76) 
(77) 



.al3 



+ 2(9"V^ + ff°V'^ 



^^J.l' 



r,ali 



q^p^ 



g^ p^ + g'^^p ]q -(g q^ + g^qjp"^, 



(78) 



where u"'^{p, q) has been defined in Eq. (j2ip . The form factors $j(s, si, S2,m?) have as entry variables, 
beside s = {p + q)'^, the virtualities of the two gluons si = p^ and S2 = q^- 

In the on-shell case only 3 invariant amplitudes contribute, which for the quark loop amplitude are 
given by 



$ig(s, 0, 0, m^ 



$2g(s, 0, 0, m2) 



$3g(s, 0, 0, n? 



9 



-.2 



+ 



g 



g'^m? 
g'^m? 



g"^ m? 
67r^s 



Co(s,0,0,m2 



2m? 



g"^ m? 
288^ ~ 247r2s2 ~ 87r2s2 



P(s,0,0,m^ 



g^ m? 
127r2s 



Co(s, 0,0, m^ 
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III 

2887r2 ' 87r2s 



+ 



g'^rin? 



+ 92Co(s,0,0,m2 



1 m 

2 + - 



4 2 n 

m m 

Att'^s 87r2 



+ ^^, ^ P(g,0,0,m^) + 



24^2^0 [s,m 



(79) 



(80) 



(81) 



247r2s 

where the on-shell scalar integrals P(s, 0, 0, m^), Co(s,0,0, m2) and Bq {s,m?) are computed in Ap- 
pendix [BJ In the massless limit the amplitude Tq'^" [p, q) takes a simpler expression and the previous 
form factors become 

(82) 



^lg(s,0,0,0) 

$2q(s, 0,0,0) 
$3q(s, 0,0,0) 



9 



367r2s' 
2887r2 s ' 



9 



2887r2 



[12L, - 35] , 



(83) 
(84) 
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Figure 5: Higher order contributions to the anomaly pole involved in the covariantization of the graviton/2- 
gluons amplitude. 



where 



L, = log( 2) *<0- (85) 

In the gluon sector the computation of Yg'^" {p, q) is performed analogously by using dimensional 



regularization with modified minimal subtraction {MS) and we obtain for on-shell gluons 

3 



rr"^(p, ?) = E ^^^(^' 0' 0) ^"' '^r"^(P' 1) 



where the form factors obtained from the explicit computation are 

11 5' 



$1^(5,0,0) 

^25(5,0,0) 



727r2s 
9' 



2887r2 s 
$39(5,0,0) = -q^Ca 



Ca, 



Ca, 



(86) 



(87) 



_^ + JLr^/S( 0) _ A^rMS(o 0) + ^ Co{s, 0, 0, 0) 



• (89) 



The renormalized scalar integrals can be found in Appendix [Bl 

The full on-shell vertex, which is the sum of the quark and pure gauge contributions, can be 
decomposed by using the same three tensor structures cf)'^^'^'^ appearing in the expansion of Tq^ {p, q) 
and T'^'"^^{p,q) 



V^^-^^^ip, q) = rf "/^(p, q) + rf -^(p, q) = Y, '^^is, 0, 0) 6'^' <^r^{p, q) , (90) 

with form factors defined as 



i=l 



$,(s,0,0) = $,,g(s,0,0) + 5^$,,g(s,0,0,m2; 



(91) 
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where the sum runs over the n/ quark flavors. In particular we have 
^2 



^i(s,0,0) 
^2(s,0,0) 



g 



727r2s 

o2 



9 ^/ z' 



Ami 



(92) 



g 



288.2, (-/-^a) 



5 



247r2 

2 



/ z' 

X^ "i? I 4 + 4^(^' 0' 0' "^?) + -^o(s, 0, 0, 



$3(s,0,0) 



5 



;iln/ - 65Ca) 

2 "/ 



1 + 



2mj 



(93) 



5^ 

87r2 



2887r2 



^e5"'(«,o) - e„"«(o,o) + sCo(s, 0,0,0) 



m,- 



- + ^V{s, 0, 0, m2) + Co(s, 0, 0, m^ 



1 + 



2mi 



(94) 



with Ca = -^'^a and the scalar integrals defined in Appendix |Bl Notice the appearance in the total 
amplitude of the 1/s pole in $1, which is present both in the quark and in the gluon sectors, and 
which saturates the contribution to the trace anomaly in the massless limit. In this case the entire 
trace anomaly is just proportional to this component, which becomes 



g 



$i(s,0,0) = -;^^(2n/-llCA) 



(95) 



The correlator T^'^""{p, q), computed using dimensional regularization, satisfies all the Ward iden- 
tities defined in the previous sections. Notice that the two-derivatives Ward identity introduced in 
Eq. 



Pc.q^T^^''''^{p,q)=0, 



(96) 



derived from the BRST symmetry of the QCD Lagrangian, is straightforwardly satisfied by the on- 
shell amplitude. This is easily seen from the tensor decomposition introduced in Eq. ([75]) because all 
the tensors fulfill the condition 



lMual3 



Paqp4>i {p,q) = 0- 



(97) 



Furthermore, we have checked at one-loop order the validity of the single derivative Ward identity given 
in Eq. (j39p and describing the conservation of the energy-momentum tensor. Using the transversality 
of the two-point gluon function Eq. (j39p this gives 



A:^r^-"^(p, 



q" p'' pl^ - q" g""^ p^ + g""/^ q'^ p"^ - g"^ p'' p ■ q) U{p 



+ [P' 



fg'^Pq^ + g'^-pPq^-g--qPp.q]li{q^), 



(98) 
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where the renormahzed gluon self energies are defined as 

n(/) = 4^fl5^P(^^o)-2 



144 7r2 



72 7r2p' 






The QCD /3 function can be related to the residue of the pole and can be easily computed starting 
from the amplitude T^'^°''^{p^q) for on-shell external lines and in the conformal limit 

<7;..r'^^"'^(p,g) = 3s$i(s;0,0,0)7x-'5(p,g) = -2^n-''(p,g), (100) 

9 



with the QCD /3 function given by 



ii?("y^^ + 3 



Pi9) = -ft2[--CA + -nf]. (101) 



As we have already mentioned, after contracting the metric tensor g^j^ with the whole amplitude F, 
only the tensor structure (f^'i"^ {p^q) contributes to the anomaly, being the remaining ones traceless, 



with a contribution entirely given by $i|m=o in Eq. (I92p . i.e. Eq. (I95p . In the massive fermion 
case, the anomalous contribution are corrected by terms proportional to the fermion mass m and 
represent an explicit breaking of scale invariance. From a direct computation we can also extract 
quite straightforwardly the effective action, which is given by 



^11CA + Inf\ j d^xd\R^^\x)U-\x,y)F^pF^P (102) 



1/3 

1 g' 

3 167r2 

and is in agreement with Eq. Q, derived from the nonlocal gravitational action. Here W^' denotes 
the linearized expression of the Ricci scalar 

R^x^ = d^^ dl h'''' -nh, h = rjf,^ h^"" (103) 

and the constant c is related to the non-abelian f3 function as 

c=-2^. (104) 

9 

Notice that the contribution coming from TJJ generates the abelian part of the non-abelian field 
strength, while extra contributions (proportional to extra factors of g and g'^) are expected from the 
TJJ J and TJJJJ diagrams (see Fig. [5]). This situation is analogous to that of the gauge anomaly, 
where one needs to render gauge covariant the anomalous amplitude given by the triangle diagram. 
In that case the gauge covariant expression is obtained by adding to the AW vertex also the AVVV 
and AVVVV diagrams, with 3 and 4 external gauge lines, respectively. 
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6 Comments 

The appearance of massless degrees of freedom in the effective action describing the couphng of gravity 
to the gauge fields is rather intriguing, and is an aspect that will require further analysis. 

The nonlocal structure of the action that contributes to the trace anomaly, which is entirely 
reproduced, within the local description, by two auxiliary scalar fields, seems to indicate that the 
effective dynamics of the coupling between gravity and matter might be controlled, at least in part, by 
these degrees of freedom. As we have just mentioned, however, this point requires a dedicated study 
and for this specific reason our conclusions remain open ended. 

Our computation, however, being general, allows also the identification of other massless contri- 
butions to the effective action which are surely bound to play a role in the physical S-matrix. They 
appear in form factors such as ^2 (Eq. [93|) and <I>3 (Eq. [M|) which do not contribute to the trace, but 
are nevertheless part of the 1-loop effective action mediated by the triangle graph. 

There are also some other comments, at this point, which are in order. Notice that while the 
isolation of the pole in the fermion sector indeed requires a massless fermion limit, as obvious from 
the structure of Eg, the other gauge invariant sector, described by Tg, is obviously not affected by 
this limit, being the corresponding form factors mass independent. This obviously does not imply 
necessarily that the gluon pole, which survives the extrapolation to the massless limit, is coupled in 
the physical S-matrix. 

Building on considerations of this nature, in particular on the possible significance of massless 
effective degrees of freedom, the role of the trace anomaly in establishing the effective interaction of 
gravity with matter has been reconsidered |101lllj. The explicit goal of this approach has been to trying 
to bypass the existing hierarchy problem between the value of the expected vacuum energy density 
{p ~ (10~^eV)^), well-described by a cosmological constant, and the Planck mass (p ~ Mp), which is 
a fundamental issue in contemporary cosmology that has not found yet a convincing explanation. In 
fact, it has been known for a long time that free massless particles contribute to the anomaly by an 
insignificant amount (T^' ~ -f^o)' proportional to the fourth power of the current Hubble rate, which 
is far too small as a value to solve the dark energy problem, due to the fact that we are living in a 
flat universe. However, it has been suggested that this small value for the vacuum energy density, 
originally attributed to the anomaly, could be raised to the expected one if the gravitational effective 
action is characterized by some effective nonlocality. In this case the contribution due to the trace 
anomaly could be modified as |19] 



Tl^^HoA^QCD-iW-'eVr, (105) 

where Aqcd is the QCD scale, which is tantalizingly close to the estimated value. While this proposal 
and similar others are clearly not the only possible solutions of the dark energy problem (similar 
values of the vacuum energy can be obtained, for instance, using axions misaligned at the electroweak 
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scale |20) and in several other ways) they share the positive feature of being characterized by few 
minimal assumptions. If so, one could envision a solution of the problem of the origin of dark energy 
without the need to enlarge the Standard Model spectrum with yet unknown particles and symmetries. 
Crucial, in these types of approaches, appears to be the role played by the effective scalar fields in 
the anomalous effective action, which are present in the local formulation of Riegert's action, together 
with their possible boundary conditions. 

7 Conclusions 

One of the standing issues of the anomalous effective action describing the interaction of a non- 
abelian theory to gravity is a test of its consistency with the standard perturbative approach. Thus, 
variational solutions of the effective action controlled by the trace anomaly should be reproduced by 
the perturbative expansion. Building on previous analysis in QED, here we have shown that also 
in the non-abelian case there is a perfect match between the two approaches. This implies that the 
interaction of gravity with a non-abelian gauge theory, mediated by the trace anomaly, indeed can be 
reformulated in terms of auxiliary scalar degrees of freedom, in analogy to the abelian case. We have 
proven this result by an explicit computation. Our findings indicate that this feature is typical of each 
gauge invariant subsector of the non-abelian TJJ amplitude, a result which is likely to hold also for 
singularities of higher order. These are expected to be present in correlators with a larger number of 
energy momentum insertions. We hope to return with a more detailed discussion of the role of the 
massless singularities - which have been found both in the trace and in the traceless part - of the TJJ 
vertex in the near future. 
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A Appendix. Feynman rules 

The Feynman rules used throughout the paper are collected here 
• Graviton - fermion - fermion vertex 
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K ( 1 



^9f.ub\h + k2)x - 2m] 



Graviton - gluon - gluon vertex 




K 

2 



-l-,5al,V^,%{kiM) 



K ( 1 



• Graviton - ghost - ghost vertex 



k2^ c^ 



(106) 



(107) 



hpu / 



* 77 ^iJ.upa k.1 p k2(T 



ki c" 
• Graviton - fermion - fermion - gauge boson vertex 

9^a f 




ig I T" Wl^^ = i^ ^ T" I - -(7^ Qua + Iv gf,a) + Qpu la 



Graviton - gluon - gluon - gluon vertex 



fl08) 



(109) 
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-g'lr'Xu'.iMkiMM) 






Graviton - ghost - ghost - gauge boson vertex 




n y J '^fjLupcr 1^2 



(110) 



fill) 



'^pupa — 9l-ip 9va I 9pu 9up 9pv 9p(J 



Dpupaih, k2) = 9pv klak2p 



9'^''kikf^ + g^p kicrk2u - gpa kip.k2u + (p ^ i^) 



Epupa{kl,k2) = 9pv {klpki^ + k2pk2a + klpk2a) 



(112) 
(113) 



gucrklpkip + g^pk2p.k2a + ifJ' ^ J^) , 

(114) 
Fpupax{kl,k2,k3) = gppga\{k2 - k3),y+gpa9px{k3 " ki)^ + g^x9pu{kl - k2)y + {fJ, o ly) 

(115) 



B Appendix. Scalar integrals 

We collect in this appendix all the scalar integrals involved in this computation. To set all our 
conventions, we start with the definition of the one-point function, or massive tadpole Ao{m'^), the 
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massive bubble BoIsjTti^) and the massive three-point function Co(s, si, S2,m^ 



Boik^m') 

Co{s, 81,82,771^) 



L^j'^'i^L^--' 


"1 

- + 1 - log 
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' /"-; ^ 


iTi^ J (P - m2) ((/ - fe)2 - m2) 


- + 2 - log ^ - 03 log , 


1 /• ,„, 


1 



7? 



(116) 



(117) 



d"/ 



ivr^ J (P — m?) ((/ — g)2 — m2"j ^^^ _j_ ^^2 _ ^^2^ 



-— y 



1=1 



Ll2 ; L12 ; h L12 ; Ll2- 



' tti + hi 



a,: - hi 



a,: - bi 



'ai + hi 



[im 



with 



1 



Am? 



8i ~v Sj -r 8}^ 



a 



(119) 



where S3 = s and in the last equation i = 1, 2, 3 and j,k 7^ i. 

The one-point and two-point functions written before in n = 4 — 2 e are divergent in dimensional 

regularization with the singular parts given by 



1 



Aoimy'^'s- ^_^2^ 



Bo{s,m 



2 \ sing. 



with 



1 1 



7 — Invr 



e e 
We use two finite combinations of scalar functions given by 

03 + 1" 



Boi8,m')m'-Aoim') = m' 



I -as log ■ 



03-1 

T>i = T>i{s, Si, m^) = Bo{s, m^) - Boisi,m'^) = 



(120) 



fl211 



fl22) 



a* log 



aj + l 
ai - 1 



03 log 



Q3 + I 
as - 1 



i = 1,2. 



(123) 



The scalar integrals Co(s,0,0, m^) and D(s,0, 0,m^) are the {si — t- 0, S2 — ^ 0} limits of the generic 
functions Co(s, si, 82, m?) and Vi{s, 81, m?) 



Co(s,0,0, m 



1 1 2 03 + 1 
■log 



Zs 03 — 1 

P(s,0,0,m2) = Vi{s,d,w?) = V2{s,d,m^) 



2-03 log 



Q3 + I 
03-1 



(124) 
(125) 
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The singularities in 1/e and the dependence on the renormahzation scale fi cancel out when taking into 
account the difference of two functions Bo, so that the Pj's are well-defined; the three-point master 
integral is convergent. 



The renormalized scalar integrals in the modified minimal subtraction scheme named MS are 
defined as 

Sp(s,0) = 2-Ls, (126) 

Sp(0,0) = -, (127) 

CO 



Co{s, 0,0,0) = - 



.^ + u;^^ + 2^^ - 


7T^^ 


~ 12. 



(128) 



where 



L, = logf-4j s<0. (129) 

We have set the space-time dimensions to n = 4 + 2w with a; > 0. The l/o; and l/o;^ singularities in 
Eqs. ()126p and ()128p are infrared divergencies due to the zero mass of the gluons. 
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